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<N ■ Abstract 

> 

(^ , The light CP even Higgs boson mass, Mh, is calculated to three- loop accuracy 

t^^ ' within the Minimal Supersymmetric Standard Model (MSSM). The result is ex- 

pressed in terms of DR parameters and implemented in the computer program H3m. 
The calculation is based on the proper approximations and their combination in 



in 

o 

f~^ ' various regions of the parameter space. The three-loop effects to M^ are typically 

of the order of a few hundred MeV and opposite in sign to the two-loop corrections. 
The remaining theory uncertainty due to higher order perturbative corrections is 
estimated to be less than 1 GeV. 
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1 Introduction 

Among the main expectations in view of the CERN Large Hadron Colhder (LHC) is to 
provide clear phenomena beyond the Standard Model (SM) of particle physics. A very 
promising candidate for an extension of the SM is the so-called Minimal Supersymmetric 
Standard Model (MSSM) which relies on an extended symmetry between fermions and 
bosons [Ill2]- It is constructed in such a way that in the low-energy limit the SM is 
recovered thus leading to the same phenomena in the energy range around the electroweak 
scale. In particular, the MSSM is in accordance with the electroweak precision data [3]. 
At the same time it provides a dark matter candidate, solves the hierarchy problem and 
provides a platform where also gravitational interactions can be included. 

An appealing feature of the MSSM is the quite restrictive Higgs sector which is described 
at leading order by two independent parameters. In particular, the mass of the lightest 
CP-even Higgs boson, M^, is not a free parameter, like in the SM, but a prediction which 
can be used in order to test this minimal supersymmetric extension of the SM. 

Mh is very sensitive to radiative corrections. In lowest order it is bound from above 
by the Z boson mass which is already excluded by experiment. Already quite some 
time ago it has been observed that large one-loop corrections, in particular from the top 
quark and top squark sector can raise Mh to about 140 GeV [4H6]. In the meantime 
a number of higher order corrections have been computed including even CP-violating 
couplings and improvements from renormalization group considerations (see Refs. [7H9] 
for a review). In this paper we consider neither CP violation nor the resummation of 
higher order logarithms. Let us nevertheless mention that in particular CP violating 
phases can lead to a shift of a few GeV in Mh, see, e.g., Refs. [?,?]. In Ref. [10] a large 
class of two-loop corrections to the lightest Higgs boson mass have been considered and 
in Ref. [H] leading logarithmic corrections at three-loop order have been computed. The 
first complete three-loop calculation of the leading quartic top quark mass terms within 
supersymmetric QCD (SQCD; more precisely, this means supersymmetric six-flavor QCD 
coupled to the MSSM Higgs sector) has been performed in Ref. [12] for a degenerate 
supersymmetric mass spectrum. It is the aim of this paper to provide details and extend 
this calculation. 

At the moment there are two computer programs publicly available which include most of 
the higher order corrections. FeynHiggs has been available already since 1998 [T3 | fT ^ [9] 
and has been continuously improved since then [151 [16]. I^i particular, it contains all 
numerically important two-loop corrections and accepts both real and complex MSSM 
input parameters. The second program, CPSuperH [TTJIIH], is based on a renormalization 
group improved diagrammatic calculation and allows for explicit CP violation. Both 
programs compute the mass spectrum as well as the decay width of the neutral and 
charged Higgs bosons. 

In this paper we discuss the three-loop corrections originating from the strong sector of 
the MSSM which are proportional to the quartic top quark mass. At three-loop order 



several mass scales enter the Feynman diagrams making their evaluation quite involved. 
In addition to the top quark mass there are the top squark masses, the gluino mass and 
the masses of the remaining squarks. An exact evaluation of the three-loop integrals is 
currently out of range. However, it is possible to apply expansion techniques for various 
limits which allow to cover a large part of the supersymmetric (SUSY) parameter space. In 
particular, we can construct precise approximations for the Snowmass Points and Slopes 
(SPS) [T9l|20]. Our set-up is easily extendable to other regions of parameter space which 
may become interesting in future. 

Together with this paper we provide a Mathematica program, HSm [21], which contains 
all our three-loop results. Furthermore, H3m constitutes an interface to FeynHiggs [22] 
and various SUSY spectrum generators which allows for precise predictions of M^ on the 
basis of realistic SUSY scenarios. 

The remainder of the paper is organized as follows: In the next Section we revisit the 
two-loop corrections. In particular, we construct approximations which are also available 
at three-loop order and compare with the exact result. In Section [3] we provide details 
on our three-loop calculations for the various hierarchies. In particular we discuss the 
renormalization and the asymptotic expansion. Section |4] describes the implementation 
of our results in the computer program H3m and the phenomenological implementations 
are discussed in Section O We present a summary and the conclusions in Section El In 
the Appendices additional material is provided, in particular all the one- and two-loop 
counterterms that have entered in our calculation. 



2 Mh in the MSSM 

2.1 Higgs boson sector of the MSSM 

The mass of the Higgs boson is obtained from the quadratic terms in the corresponding 
potential which for the MSSM has the following form: 

Vh = (I/UsusyI^ + mf) \Hif + (I/xsusyI^ + ml) \H2\'^ 

— "^12 (^abHiH2 + eabH^*H2 



+ l{9l + 9l)[m'-m'Y + \gl\HlH^. 
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with ei2 = —621 = 1 and en = €22 = 0. /xsusy is the Higgs-Higgsino bilinear coupling 
from the super potential and mi, 7712 and rriu are soft breaking parameters. Note that the 
quartic terms are fixed by the SU(2) and U(l) gauge couplings gi and (72- The parameters 



in (11]) are related to the masses of the gauge bosons and the pseudoscalar Higgs via 

Ml = -ml^ (tan/3 + cot (3) , 

where tan/3 = f2/fi. 

After spontaneous symmetry breaking the neutral components of the doublets Hi and H2 
acquire the vacuum expectation values vi and V2, and we write 

Hi= A'l + TI^'^^ + '^i' 
V -01 

which leads to the following representation of the Higgs boson mass matrix at tree level: 

2 _ 1 d'^V _ sin 2/3 / M| cot /3 + Ml tan /3 -M^ - M^ \ 

-^H.trco - 2 50^502 ~ 2 V -Ml -Ml M|tan/3 + Mlcot/3 J '^ ^ 

where we restrict ourselves to the two CP-even Higgs bosons. 

In this paper we restrict ourselves to the numerically dominant mf corrections. As a 
consequence we can set the electroweak gauge couplings to zero and furthermore nullify 
the external momentum in the occurring two-point functions. The formalism presented 
in the following is adapted to this framework. 

It is convenient to evaluate the quantum corrections in the {0i,02} basis which requires 
the evaluation of self energy corrections E^^ and E<^2 involving (pi and 02. Denoting 
renormalized quantities with a hat, Eq. (jll) gets modified to 

Mj, = Ml,,,, - ( ^'^ ^p^^ ) , (5) 



with [15 



T.^^ = S<^j - Sa sin^ P 



+ ^TJJ — ^^^</>i cos /3 (1 + sin^ /?) 



e ■ .2 



t02 COS /3 sin /3 . 



2MvKsin^?vF ' 
E<^2 = E^j - Ea cos^ /3 

" ^m — ^^^01 sin^ /3 cos /3 
+ ^777 — ^^—^02 sin /3 (1 + cos^ /3) , 
]</,i^2 = E<^^<^2 + Ea sin /3 cos /3 (6) 

-t^i sin^ /3 
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In this equation, ■i^vi/ is the weak mixing angle, E^ denotes the self energy of the pseudo- 
scalar Higgs boson and t^. the tadpole contributions of the field 0j. Typical diagrams to 
the individual contributions can be found in Fig. [TJ Since we are only interested in the 
leading corrections proportional to m^ we evaluate the quantum corrections in the limit 
of vanishing external momentum [231425] . 

For the evaluation of the lightest Higgs boson mass we consider in a first step the matrix 
element of A^|^ to a given order in perturbation theory. Subsequently, we determine the 
eigenvalues and assign the smaller one to M^. We perform this procedure at tree level 
and at one-, two-, and three-loop order which leads to the corresponding approximations 
of the Higgs boson mass. It is convenient to introduce the quantity 

AMf = M^'-'°°P) - M*'^^% (7) 

representing the difference between the Higgs boson mass evaluated with i-loop accuracy 
and the tree-level result. 



2.2 Top squark sector of the MSSM 

In this paper we only consider strong corrections which means that apart from the quarks 
and gluons also the corresponding superpartners, the squarks and gluinos, are present. 
The leading contribution proportional to Gpm^ is generated by the top quark Yukawa 
coupling which distinguishes the top squark sector from the other squark parts of the 
MSSM Lagrange density. In order to fix the notation let us discuss in more detail the 
mass matrix of the left- and right-handed component of the top squark, tL and in, which 
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Figure 1: Sample diagrams contributing to S^^, T,^^, ^<i>i<t>2^ ^A, i(t>i and t^^. Internal 
solid, dashed, dotted and curly lines correspond to top quarks, top squarks, 5-scalar and 
gluons, respectively. Gluinos are depicted with as curly lines with an additional solid line 
in the middle. The external dashed line corresponds to the Higgs boson. 



is given by 

_ /m? + M| (I - I sin^ dw) cos 2/3 + M? 
* V ^t {At - /isusY cot (3) 
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TUt {At - /xsusY cot 13) 



ml + lMl sin^ ^w cos 2/3 + M? 



with Xt = At — /isusY cot /3. Mq and M^ are soft SUSY breaking masses, and At is the 
soft SUSY breaking tri-linear coupling between the Higgs boson and the top squark fields. 



Diagonalization of Eq. (|8]) leads to the mass eigenstates ti and ^2 with masses 
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The mixing angle is defined through the unitary transformation 

"j^ \)=n\Mint. with 7^,= f^°^^* "^^"/^ 

m~ J t t ti ' \^ sm ^t cos b^t 



(10) 



and 
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2.3 Leading m^ corrections in the on-shell and DR scheme 

As already mentioned above, the numerically dominant contribution arises from the self 
energy diagrams evaluated for vanishing external momentum. Thus, it is convenient to 
introduce the following notation for the z-loop corrections to the Higgs boson mass 

AM^ = A'^'Mf +A''"'°Mf , (12) 

where A"^*M^ comprises the complete SQCD contribution of order ata^~^ {at is the top 
Yukawa coupling) originating from the top quark/squark sector for vanishing external 
momentum which is proportional to m^. At one-loop order only top quarks and top 
squarks are present in the loops. The two- and three-loop corrections. A™* M^ and 
A"** M^ , are obtained by adding gluon, gluino, quark and squark contributions. A'''^™M^ 
represents the remaining part which is only available at one- and two-loop order. In 
our approach these corrections are taken from FeynHiggs which includes the complete 
one-loop corrections and all available two-loop termsjj At three-loop order only the 
contribution A"^* M^ is considered. 

In the following we discuss the relative contribution to the Higgs boson mass comparing 
AM^ and A^'tM^ at one- and two- loop order {i = 1,2) where both the on-shell and 
DR scheme for the mass parameters and the mixing angle are considered. For illustration 
we adopt the scenarios SPSla and SPS2 and show the Higgs boson mass as a function of 

mi/2. 

The solid lines in Fig. [2] show the results for AM^ using on-shell parameterqj for the 
masses and 6t as it is provided by FeynHiggs. The dashed lines correspond to the mf 
approximation (A'^tM^ ) which can be found in Refs. [T51I26] and has been confirmed by 
us by an independent calculation. The panels (a) and (b) correspond to the SPSla and 
SPS2 benchmark scenarios, respectively, where the input parameters have been generated 
with the help of SOFTSUSY [27|. The small differences of the solid and dashed linecl 
demonstrate that the leading top quark mass term approximates the full result to a 
high accuracy. This statement is also true in the DR scheme. The corresponding two-loop 
results are shown in Fig. [2] as dash-dotted and dotted curves where the former corresponds 
to the full and the latter to the approximate result. 

It is well known that the perturbative series can exhibit a bad convergence behaviour in 



^For a detailed description we refer to the FeynHiggs home page [55]. 

^For the on-shell renormalization of the mixing angle we adopt the convention of Ref. 

■^In the case of SPSla the two-loop dashed line is almost on top of the solid one. 
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Figure 2: Comparison of complete and approximate one- and two-loop corrections to the 
Higgs boson mass for SPSla (a) and SPS2 (b). The solid (full result) and dashed lines {mf 
approximation) represent the results in the on-shell scheme where the upper and lower 
curves correspond to the one- and two-loop results, respectively. The two-loop DR results 
are shown as dash-dotted (full result) and dotted (m^ approximation) curves. 



case it is parametrized in terms of the on-shell quark masseqj which is due to intrinsically 
large contributions related to the infra-red behaviour of the theory. Thus, it is tempting 
to re-parametrize the results for the Higgs boson mass in terms of DR parameters for 
the top quark mass, the masses of the SUSY particles and the top squark mixing angle. 
Since the corrections are dominated by A"** M^ it is sufficient to consider only this term 
and take over A'^'^^M^ from the output of FeynHiggs. In the remainder of the paper we 
will refer to this renormalization scheme as DR scheme although it contains a mixture of 
on-shell and DR parameters. Let us mention already at this point that from outside this 
mixture does not pose any complication in the practical use since the spectrum generator 
produces both on-shell and DR parameters which then serve as input for the evaluation 

Further below we will discuss the scheme dependence of Mh and indeed show that the 
loop corrections are in general smaller in the DR scheme (cf. Fig. [8]). Similar studies can 
also be found in the literature [T^IH]. 

The considerations of this subsection motivates the following procedure at three loops: It 
is certainly sufficient to consider only the approximation A™* M^ since the size of the 
remaining term is expected to be below 100 MeV. Furthermore, we adopt the DR scheme 
since we expect that the perturbative series shows a better convergence behaviour. In 
addition the evaluation of the counterterms themselves is significantly simpler. Actually, 
most of them are already available in the literature and the computation of the remaining 
ones is quite straightforward as we discuss in Appendix A. We provide the analytical 



"^For a typical example we refer to the electroweak p parameter. Using the on-shell top quark mass 
the four-loop corrections [28H30] are larger by a factor 50 as compared to the MS scheme. 



results for the two-loop renormalization constants of the top squark masses and mixing 
angle, that can be easily expanded for the mass hierarchies considered in this paper. The 
multiplicative DR renormalization constants of the top quark and gluino mass are mass 
independent and therefore valid for all hierarchies. 



2.4 Construction of approximations 

Considering the many different mass parameters entering the formula for the Higgs boson 
mass an exact calculation of the three-loop corrections is currently not feasible. However, 
due to the various hierarchies among the particle masses it is promising to consider ex- 
pansions in properly chosen small parameters. As a guideline for the latter we follow the 
SPS scenarios as defined in Refs. [T91I20]. 



In order to construct approximations covering all SPS cases it is sufficient to consider the 
following hierarchies among the SUSY masseCl 



(h3) 


rrig ^ 171^^ ^ mi^ ^ m-g , 


(h4) 


mq > mi^ K. mi^ ^ m-g , 


(h5) 


niq > mi^ > mi^ k. m-g , 


(h6) 


mg > mi^ ^ m-g:^ m^^ , 


(h6b) 


mg ^ m^^ ^ mg :^ m^^ , 


(h9) 


mg ^ m^^ ^ mi^ > m-g , 



(13) 

where in the case of "S>" an asymptotic expansion in the corresponding hierarchy is 
performed. In the case of "~" a naive Taylor expansion in the difference of the particle 
masses is sufficient. Throughout this paper, q denotes any squark other than t, and we 
assume a common mass value nig = m,g-^ = m,q^ for all of these "heavy squarks" . 

In all hierarchies we assume that the SUSY masses are larger than the top quark mass 
and perform an asymptotic expansion in the corresponding ratio. In the numerical results 
discussed below we include for the various hierarchies the expansion terms as given in 



■'We decided to keep the nomenclature for the hierarchies as they are in our internal computations 
and documents. The non-continuous numeration results from the fact that for testing purposes we have 
computed further hicrachics which, however, are not included in the program H3m (cf. Section |4l). 
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Table 1: Expansion terms available for the individual hierarchies as defined in Eq. ( I13p 
at three-loop order. 



Tab. [T] where the following notation has been introduced 
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Note that at two-loop order the contributions involving the squarks q with q G {u, d, s, c, b} 
cancel in the sum of all diagrams. At three-loop level, however, the results depend on rriq. 
In those cases where rriq is much larger than the other masses at least three expansion 
terms are computed and a good convergence even up to rriq ^ ttt-j^ is observed. 

In order to demonstrate this point we consider the hierarchies (h3) and (h4) and show in 
Fig. [3] the three-loop prediction for Mh. The dashed and solid lines correspond to (h4) 
including successively higher orders in l/ruq where for illustration the following input 
parameters have been chosen|j 
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The horizontal (??7,g- independent) dotted line corresponds to the scenario (h3) with rriq = 
^susY fixed at 800 GeV. One observes a crossing of the latter and the (h4)-curve including 
1/m^ corrections for rrig ^ t^susy which nicely demonstrates the rapid convergence in the 



^If not stated otherwise we set the renormahzation scale equal to the on-shell top quark mass and 
evaluate all DR parameters at that scale. Note, however, that H3m is not restricted to this choice. 
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Figure 3: Dependence of the three-loop corrections on the heavy squark mass rriq. The 
dashed and solid lines correspond to the hierarchy (h4) where successively higher or- 
der terms in msusY/f^q have been included (The solid curve contains terms of order 
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The dotted and dash-dotted curves correspond to the hierarchy (h3); for 



the dotted line rriq has been kept fixed at 800 GeV, the dash-dotted curve includes terms 
of order (1 — 'm'^^^y / m'^Y . 



large- 771^ expansion. Fig. 12] also shows that the expansion around rriq = rn^vsY leads to 
good approximations even if rriq is two to three times as big as msvsY, see dash-dotted 
curve. Note that for this plot we computed the DR top quark mass for rriq = 800 GeV 
and kept it fixed. 

Let us in a next step compare the approximate SQCD corrections according to our hi- 
erarchies with the full prediction for Mh from Ref. [26]. For illustration we adopt in the 
remainder of this Section a minimal supergravity (msugra) scenario with 
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Figure 4: Comparison of approximate and full two-loop result for hierarchy (h3) (a) and 
the combination of (h3), (h5), (h6), (h6b) and (h9) (b). The contour lines indicate the 
deviations in MeV. In (c) and (d) the results of (a) and (b) are normalized to the genuine 
two- loop contributions where the contour lines indicate the deviations in per cent. The 
benchmark points and slopes are shown as (white) dots and lines. 



and vary mo and mi/2 as follows 



60 GeV <mo< 1600 GeV, 
100 GeV < mi/2 < 800 GeV . 



(17) 



We have checked that very similar results are obtained for other choices of tan/3, Aq, and 
sign(/isusY)- Thus, our conclusions are at least valid for all MSUGRA SPS scenarios (cf. 
Appendix C). 

In Fig. Ht^a) and (b) the absolute value of the difference between the full and the ap- 
is shown in 
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proximate two-loop prediction for the Higgs boson mass, M^ 
the mQ-mi/2 plane where M^ ™ includes the same number of expansion terms which 
are available at three loops. In Fig. Hl^a) we only include the results from hierarchy (h3) 
whereas in (b) also (h5), (h6), (h6b) and (h9) enter. For each hierarchy we compute the 
difference to the exact result and plot in Fig. Hl^b) the minimum. We define the relative 
uncertainty through 
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Table 2: Comparison of full and approximate two-loop prediction for M^ for the different 
benchmark points. 



which is shown in Fig. HJ^c) and (d). In Eq. f ll8p M^ corresponds to the exact z-loop 
prediction. For reference we show in Fig. |4] the MSUGRA SPS benchmark points and 
slopes as (white) dots and lines, having in mind that for some of them the values of tan /3 
and Aq are different from the ones chosen in Eq. (TT6|) . The assignment of the individual 
scenarios to the corresponding dot is easily done with the help of the table in Appendix C 

Already for (h3) alone one observes a good coverage in the whole mQ-mij2 plane with 
deviations smaller than 150 MeV. This gets further improved after including the other 
hierarchies. For lower values of mi/2 one has relative deviations also above 10%, however, 
the absolute difference between the full result and the approximation is below 100 MeV. 

In Tab. [2] we directly compare the two-loop predictions for Mh for the benchmark points 
listed in the table in Appendix C and SPS7 and SPS8 (gauge-mediated supersymmetry 
breaking) |T9]. As before, the full results are based on FeynHiggs and Ref. [26], and 



for M^ -''^PP'' ^g ygg ]^}^Q approximation incorporated in HSm. An impressive agreement is 
found, often even below 100 MeV. 

The results discussed in this Subsection are very promising in view of the three-loop 
approximation. At two-loop order the expansion terms specified in Tab. [1] provide an 
excellent approximation to the full result. Thus, it can be assumed that the corresponding 
terms at three loops approximate the unknown result with high precision. 



3 Technical details to the three-loop calculation 

The three-loop calculation of the individual Green's functions contributing to Mh is or- 
ganized as follows: All Feynman diagrams are generated with QGRAF [31]. In order to 
properly take into account the Majorana character of the gluino, the output is subse- 



ts 



quently manipulated by a PERL script [32] which apphes the rules given in Ref. [33]. The 
various diagram topologies are identified and transformed to FORM [3l| with the help of 
q2e and exp [35l|36]. The program exp is also used in order to apply the asymptotic 
expansion (see, e.g., Ref. [37]) in the various mass hierarchies. The actual evaluation of 
the integrals is performed with the package MATAD [38] , resulting in an expansion in ci — 4 
for each diagram, where d is the space-time dimension. 

The total number of three-loop diagrams amounts to 6706 and 7670 for the (pi and 02 
self energies, respectively, and 845 and 982 for the corresponding tadpole contributions. 
The computation of the off-diagonal matrix element S<^^2 involves 6136 diagrams, and the 
propagator of the pseudoscalar Higgs another 7670. The application of the asymptotic 
expansion significantly enlarges these numbers leading to about 37000 (for (h3)) or even 
94 000 (for (h6)) subdiagrams. Note that there are diagrams where, depending on the 
hierarchy, up to 15 subdiagrams have to be considered. A typical example is shown in 
Fig. [la). 

A subtlety arises from diagrams as the one shown in Fig. [T](b). If both the external mo- 
mentum and the e-scalar mass are set to zero from the beginning, an infra-red divergence 
occurs and cancels the ultra-violet divergence of the integral. In effect, the diagram will 
be of order {d — 4) due to the e-scalar algebra. In order to avoid this, we keep the external 
momentum q non-zero, though much smaller than all other scales. The ultra-violet pole 
multiplied by the algebraic factor of {d — 4) then produces a finite contribution, while the 
infra-red divergence leads to [d — 4) In(g^) and vanishes as d — )■ 4. 

Instead of the requirement g 7^ one could also introduce a nonzero mass for the e-scalars 
in order to regulate the infra-red divergences. In the final result we again observe that 
the regulator is multiplied by an additional factor [d — 4) leading to a finite result for 
Me — )■ 0. We have checked that the latter prescription leads to identical results as the 
one with g 7^ 0. 

We refrain from presenting all available analytical results for the expansion in the various 
regions. They are implemented in the program H3m and thus easily accessible if necessary. 
However, for the convenience of the reader we provide in this Section the result for (h4), 
see Eq. (TT3I1 . which could be useful for other applications. We will present the results 
expressed in terms of the DR parameters a^, mt-, rn^^, mi^, rrig and Ot- The corresponding 
counterterms can be found in Appendix A. 

Before providing explicit expressions a comment concerning the DR renormalization con- 
stants for the top squarks is in order. Due to diagrams involving heavy squarks q, for 
example Fig. Ej^a), the squared Higgs boson mass receives contributions which are propor- 
tional to rrig and thus can lead to unnatural large corrections. For this reason we adopt 



the on-shell scheme for these contributions to SZm^ and SZ^^ (cf. Eq. f l37|) and Fig. [5t^b) 
for a sample diagram). This avoids the potentially large terms ~ ?7i| from the three-loop 
diagrams. We follow this procedure also in the case where the top squarks and the "heavy 
squarks" are degenerate in mass. The renormalization of the mixing angle is free of such 
enhanced contributions and we can stick to the pure DR scheme in that case. We hasten 
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Figure 5: (a) Feynman diagram involving a heavy virtual squark contributing to the 
Higgs boson self energy, (b) Counterterm diagram related to the diagram in (a). The 
same notation as in Fig. [1] has been adopted. 

to add that this discussion only concerns the internal structure of H3m and has no direct 
consequences for the user. The input parameters of H3m are the DR ones as they appear, 
e.g., in the output of SOFTSUSY. 

As already noticed in Refs. [I5l|26], a similar behaviour is observed when the gluino is 
much heavier than the top squarks. In this case, the two- and three-loop corrections 
to the Higgs masses computed in the DR scheme contain terms proportional to rra^ and 
771? These contributions are canceled in the on-shell scheme by the finite parts of the 
relevant counterterms. Thus, in order to avoid unnatural large radiative corrections to the 
Higgs masses, we adopt for scenarios with heavy gluino masses a modified renormalization 
scheme for the top squark masses. We call this scheme "modified DR" (MDR) and it 
is characterized by the non minimal renormalization of the top squark masses. The 
additional finite shifts of top squark masses are chosen such that they cancel the power- 
like behaviour of the gluino contributions. Again, the renormalization of the mixing angle 
will not be modified as compared to the genuine DR scheme. 

The relevant finite shifts for the scenarios considered in this paper are explicitly given in 
Appendix B. As can be noticed from Eq. (IT3ll the scenarios (h4), (h5) and (h6) display 
heavy squark mass contributions whereas heavy gluino terms are specific only for the 
scenarios (h6) and (h6b). 

In the practical calculation we use the DR top squark mass parameters as provided by the 
spectrum generators and transform them with the help of the formulae of Appendix B 
to the corresponding parameters in the MDR scheme which constitute the input for our 
analytic expressions. 

In the following we present results for the renormalized two-point functions S^-^ , S^j and 
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^We only include terms up to order 1/m?. 
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with mi = mtinr), msusY = "^susYl/^r) = m^il^r) = m-t^M = "^gl/^r), /^t = ln(/i?/m2), 
/i5 = ln(m^/mgugY) and Itg = ln(m^/?Ti?) where fir is the renormahzation scale. The 



on-shell result corresponding to T,^^ has been presented in Ref. [12] for At = 0. 

We refrain from providing more analytic results since all of them come along with the 
program H3m which is discussed in the next section. 



4 Description of H3m 

In this Section we describe the implementation of our three-loop results in a user-friendly 
computer program which allows the evaluation of the light CP even Higgs boson mass 
Mh to three- loop accuracy. The program is implemented in the form of a Mathematica 
package. 

To set the input parameters for the calculation, i.e. the SUSY spectrum and SM param- 
eters, the SUSY Les Houches Accord (slha) [52] is used. For ease of use, we provide 
functions that call a spectrum generator from Mathematica to produce an SLHA spec- 
trum file. To produce the plots in this publication, we have chosen SOFTSUSY [27], but it 
is possible to use any spectrum generator that provides the DR parameters in addition to 
the on-shell mass spectrum like SuSpect [ID] or SPheno [H]. The advantage of SOFTSUSY 
is that the renormahzation scale of the DR parameters can be chosen independently of 
the electroweak symmetry breaking. 

The corrections A™* Mh to Mh, being proportional to the fourth power of the mass of the 
top quark rrit, are very sensitive to both the definition and the uncertainty of rrit. Thus, 
it is important to use the most precise value of rrit available. For this reason we take into 
account the full two-loop SQCD corrections between the on-shell and DR top quark mass 
given in Ref. [12] |j In this paper, the relation of the on-shell top mass Mj and DR top 
mass rrit was derived as a function of the DR masses. Solving this equation iteratively. 



^We thank Steven Martin for providing us with the relevant formulae from Ref. |42j in electronic form, 
and for allowing us to include his code in our program. 
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we get nit as a function of Mj. The integrals appearing in |12] are evaluated using the 
C library TSIL [l3]. This relation is available for general renormalization scale Hr which 
enables us to obtain mt{jJir) in the DR scheme using the on-shell mass Mt as measured at 
the Tevatron [H] as input. 

Another critical parameter for the evaluation of M^ is the strong coupling a^. We use 
as{Mz) as input and follow Ref. [ISHIT] in order to evaluate a^ in the DR scheme with all 

SUSY particles contributing to the running. First, ah (Mz) = 0.1184 [IH] is run up to 
the decoupling scale, which we set to the average value of the SUSY particles, using the 
four-loop P function [l9l[50]. There, we perform the transition to the DR scheme and the 
full theory. The two-loop matching coefficients from |17] are used in this step. To obtain 

Q!s (fXr) for arbitrary values of the renormalization scale fir, we use the three-loop 

SQCD 13 function given in |SI1I32]. 

The remaining input parameters comprise the ones for the SUSY breaking scenario, which 
we summarise in the table of Appendix C and the SM parameters Mz, Gp and a which 
also serve as input for the spectrum. The default values set in H3m read 



Mz -- 


= 91.1876 GeV, 


Mt = 


= 173.1 GeV, 


Gf -- 


= 1.16637-10"^ GeV~2 


/a{Mz) -- 


= 127.934, 


as{Mz) -- 


= af '^S(Mz) = O.m 



(19) 

Of course, it is possible to modify these default values. Note that the top squark masses 
and mixing angle are obtained from the soft breaking parameters according to Eqs. iQ 
and (ITT]). 

In order to include all the known corrections to Mh at the one- and two-loop level, the 
spectrum file is passed to FeynHiggs P| [T3HT5] . FeynHiggs uses on-shell parameters that 
are given in the spectrum file and provides the neutral Higgs mass matrix up to the 
two-loop level. Since we prefer to use the DR scheme, we need to perform a conversion 
before adding our three-loop results. This we do by subtracting the on-shell expression 
for A"** Mh (up to two loops, without any expansions in the masseqj) and adding it back 
in the DR scheme. Thus, we use the DR scheme to evaluate A"^* Mh, which are dominant 
and sensitive to the top quark mass, and the on-shell scheme for A'''^™M/j. 

The next step is to choose a suitable mass hierarchy for the expansion of the three- 
loop corrections. This is done by comparing, at the two-loop level, the full result from 
Ref. [2S] with the expansions in all the mass hierarchies and choosing the one minimizing 
the error. Finally, the three-loop corrections are added, the neutral Higgs mass matrix is 
diagonalized, and the mass of the light Higgs is returned to the user. 

The interface of the program is outlined in Fig. |6l The parameters are set up by a call of 



'We thank Pictro Slavich for sending us the compact formulae from Ref. [26] in electronic form. 
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Figure 6: Flowchart of H3m. First, the user calls HSSetSLHA or one of its descendants to 
set the parameters. A subsequent call to HSmcomputes M^. 

the function HSSetSLHA, which passes its arguments to the spectrum generator and parses 
its output to get the relevant input parameters for the calculation. The top mass and 
strong coupling are calculated as described above. Alternatively, the function HSGetSLHA 
uses an existing spectrum file instead of running a generator. For the user interested in 
the Snowmass Points and Slopes, we provide convenient wrapper functions H3SetSPS<x> 
which call HSSetSLHA with parameters according to a specific benchmark scenario. 

The main calculation is organized by the function HSm, which calls FeynHiggs, does the 
conversion to the DR scheme described above, chooses an appropriate mass hierarchy, 
adds the three- loop corrections, and returns M/j. 

In Fig. [TJ a typical JVIathematica session with HSm is shown. A more detailed description 
of HSm comes along with the program which can be found at the web page [2Tj . 



5 The Higgs boson mass to three-loop accuracy 

In this Section we use, if not stated otherwise, the input parameters as listed in Eq. ( TTOl) 
and furthermore adopt for the renormalization scale /i,- = Mt as our default value. 

In Fig. |H]we show the renormalization scheme dependence of Mh as a function of mi/2 for 
the SPS2 scenario. This is convenient since we have the same abscissa both for the on-shell 
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Mathematica 7.0 for Linux x86 (64-bit) 
Copyright 1988-2008 Wolfram Research, Inc. 

In[l] := Needs ["H3"'] ; 

RunDec: a Mathematica package for running and decoupling of the 

strong coupling and quark masses 
by K.G. Chetyrkin, J.H. Kuhn and M. Steinhauser (January 2000) 

In [2]:= H3SetSPSla[ 300.]; 

H3GetSLHA: :TSIL: Using TSIL by S.P. Martin. 

FeynHiggs 2.6.5 

built on Dec 20, 2008 

T. Hahn, S. Heinemeyer, W. Hollik, H. Rzehak, G. Weiglein 

\protect\vrule widthOpt\protect\href -[http : //www . f eynhiggs . de}-[http : //www . f eynhiggs . de} 



FHHiggsCorr contains code by: 

P. Slavich et al . (2-loop rMSSM Higgs self -energies) 
Loading Results for hierarchy h3 
Loading Results for hierarchy h3 
Loading Results for hierarchy h6b2qg2 
Loading Results for hierarchy h6b2qg2 

In [3] := H3m[] 

Loading Results for hierarchy h6b2qg2 

Out [3]= {mh -> 114.176} 

Figure 7: A typical Mathematica session with H3m. 

and DR result. Note, however, that the three- loop on-shell result is only available for a 
degenerate mass spectrum of the SUSY particles and vanishing parameter At [12]. Thus, 
we restrict ourselves to this limit also for the DR result. For this reason the following 
discussion should be considered in a less quantitative but more qualitative sense and 
should not be used, e.g., for estimating a theoretical uncertainty. 

In the left panel of Fig. |8]the upper dotted, dashed and solid curve correspond to the one-, 
two- and three-loop prediction of Mh in the on-shell scheme whereas the corresponding 
lower three curves are obtained in the DR scheme. In the on-shell scheme one observes 
large positive one-loop corrections which get reduced by 10 to 20 GeV after including the 
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Figure 8: Renorinalization scheme dependence of M^ as a function of mi/2 adopting SPS2. 
Dotted, dashed and sohd curves correspond to one-, two- and three-loop results. The DR 
(on-shell) results correspond to the lower (upper) three curves. In the right panel the 
interesting part of the left one is magnified. 

two-loop terms. The three-loop corrections amount to several hundred MeV. They are 
positive or negative — depending on the value of mi/2. 

The situation is completely different for DR mass parameters: the one-loop corrections are 
significantly smaller and lead to values of M^ which are already of the order of the two- 
and three-loop on-shell prediction. The two-loop term leads to a small shift of the order of 
— 1 GeV and the three- loop term to a positive shift of about the same order of magnitude. 
The final prediction for Mh is very close to the one obtained after incorporating three-loop 
on-shell result so 

Comparing the DR and on-shell results in Fig. [8] one observes a nice reduction of the 
scheme dependence when incorporating higher order corrections]^^! Whereas there is a 
huge gap between the two one-loop curves (dotted) the difference in the two-loop pre- 
diction of Mh is below 2 GeV which gets further reduced by about a factor ten after 
incorporating the three-loop corrections to roughly 200 MeV. 

Fig. [H] extends Fig. IHto three loops. In (a) and (b) we again discuss the hierarchy (h3) and 
the combination of the hierarchies (h3), (h5), (h6), (h6b) and (h9), respectively, and show 
the difference between our best prediction and the one where the expansion parameters 
are cut by one unito In the whole parameter plane we observe small absolute corrections 
reaching at most about 100 MeV. This leads to the conclusion that as a conservative 
estimate of the uncertainty of our approximation procedure one can assign about 100 MeV. 



^°The relatively large three-loop corrections (as compared to the two-loop ones) do not pose any 
problem since we use simplified formulae as mentioned before. Furthermore, there are regions in the 
parameter space where the two-loop corrections are accidentally small in the DR scheme leading to large 
relative three-loop terms. Nevertheless the overall size of the two- and three-loop corrections is small. 
^^Up to two- loop order the scheme dependence has already been discussed in Fig. [21 
^^ We only cut in parameters originating from asymptotic expansion: when counting powers of mass 
ratios, we leave sin 26*4 untouched and do not replace it by Eq. (|TT|) . 
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Figure 9: Comparison of the three- loop predictions for Mh using the maximal available 
expansion terms and reduced input as described in the text. In (a) and (b) the absolute 
deviations are shown for the hierarchy (h3) and the combination of (h3), (h5), (h6), (h6b) 
and (h9), respectively. The contour lines indicate the deviations in MeV. In (c) and (d) 
the results of (a) and (b) are normalized to the genuine three-loop contributions where 
the contour lines indicate the deviations in per cent. The benchmark points and slopes 
are shown as (white) dots and lines. 



In Fig. M,c) and (d) we show relative deviations defined through 



^(3) 



Mf ) - Mf ^''''* 



M, 



(3) 



(2) 



r(3) 



M, 



(3), cut 



(20) 



where M^ is our best three- loop prediction, M^ is the prediction where at three- 

loop order the expansion depth of each parameter (cf. Tab. [T]) is reduced by one unit, 
and M^ corresponds to the (full) two-loop term. Similar to the two-loop case, larger 
corrections are only observed for small values of mo and mi/2 which is a consequence of 
a small denominator in Eq. ( 120|) . The three-loop correction terms, however, are stable as 
can be seen from the panels (a) and (b). Thus, as in the two- loop case, we are able to 
cover the whole mo-mi/2 plane and are in particular able to produce precise values for 
Mfi for all SPS scenarios. 

The slopes for three SPS scenarios are shown in Fig. fTOT a). (b) and (c) where the dotted, 
dashed and solid lines correspond to the one-, two- and three-loop predictions, respectively. 
For all three cases one observes negative corrections between 1 and 4 GeV at two loops 
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Figure 10: M^ for the different slopes of tlie benchmark scenarios SPSla (a), SPS2 (b) 
and SPSS (c). Dotted, dashed and sohd hnes correspond to the one-, two- and three-loop 
predictions. The dashed line with longer dashes (at two loops) correspond to the full 
results, the one with the shorter dashes to the approximation implemented in H3m. In (d) 
the dependence of M^ on the renormalization scale is shown where the dotted, dashed 
and solid line corresponds to the one-, two- and three-loop prediction. 

and positive contributions from the three-loop term which amount up to about 2 GeV. 

In Fig. [TUT d) we show the dependence of the prediction for M^ on the renormalization scale 
/i^. As an example we adopt the SPSla scenario with mi/2 = 250 GeV and exploit that 
SOFTSUSY allows the evaluation of all DR parameters at the scale [i^. One observes a strong 
dependence at one-loop order which gets significantly reduced at two loops. The three- 
loop curve is even more flat resulting in a stable prediction for M^. Around /i^ = 250 GeV 
the two-loop correction shows a local maximum and is furthermore very small whereas 
the three- loop term still amounts to about 500 MeV. Around /i^ = ^t, which is often 
used as a default choice, one has negative two-loop corrections of about —2 GeV and a 
slightly larger three-loop contribution than for /i^ = 250 GeV. The corresponding plots 
for SPS2 and SPSS look very similar. Thus, we refrain from presenting them here; they 
can easily be generated with the help of HSm. 
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Figure 11: Dependence of Mh on At. The dotted, dashed and solid line corresponds to 
the one-, two- and three-loop prediction. 

It is interesting to investigate the dependence of Mh on the soft breaking parameter At. 
In Fig. [TT] we show the result for M^ where the following values for the parameters have 
been chosen 



TTlf 



m 



t2 



m-r, 



500 GeV , 
1000 GeV , 
500 GeV , 



rriq = 2000 GeV , 

/isusY = 800 GeV, 

tan /3 = 10 , 

Ma = 1500 GeV . 



(21) 



Furthermore, we employ only the mf corrections since it is not possible to transmit the 
parameters of Eq. ( ETj) directly to FeynHiggs and evaluate the corresponding Higgs boson 
mass. 

It is interesting to note that the three-loop corrections are quite sizeable, amounting 
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Figure 12: Genuine two- (upper panel) and three-loop (lower panel) corrections to M^ in 
the mQ-mii2 plane. 

up to about 3 GeV. In contrast to the two-loop terms they are positive and lead to a 
compensation. For At = — 2TeV and At = the three-loop prediction is even above the 
one- loop value for Mh. 

In order to get an impression on the size of the three-loop corrections we show in Fig. [12] 
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Figure 13: Prediction of M^ to three-loop accuracy using H3m. The same conventions as 
in Fig. [12] have been adopted. 

(lower panel) the difference between our best three-loop prediction and the full two-loop 
result as a function of ttlq and rrai/2. We observe that the corrections are always positive 
and vary for our parameters between a few hundred MeV and about 2 GeV. They show 
only a mild dependence on mo, but vary strongly with mi/2- In particular the corrections 
become larger for increasing values of ?7ii/2- For comparison, we show in Fig. [12] (upper 
panel) the corresponding quantity at two- loop order, i.e., the difference between the two- 
loop and the one-loop result. In contrast to the three- loop contributions they are negative 
and amount to about twice the three-loop terms in a large region of the parameter space. 
However, there are also regions where the three-loop corrections are larger than the two- 
loop ones. Note, however, that in the whole mo-mi/2 plane the one-loop corrections are 
more than ten times bigger. Furthermore, the occurrence of three-loop corrections, which 
are large compared to the two-loop ones, can also be seen from Fig. [TOT d): whereas the 
two-loop corrections vary between —4 GeV and +4 GeV in the considered range of fir the 
three-loop term is almost constant and amounts to about 1 GeV relative to the two-loop 
result. 



In Fig. [13] we finally show the three-loop prediction of Mh including the three-loop effects 
discussed in this paper. Again we restrict ourselves to the parameter space defined above. 
Values for other input parameters are easily obtained with the help of H3m. One can 
see that for increasing mi/2 also the Higgs boson mass gets larger and values well above 
120 GeV can be reached. This is already observed at one-loop order and is due to the 
fact that for larger values of mi/2 the whole supersymmetric spectrum becomes heavier. 
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Mf 


^(3),cut 


optimal 




(GeV) 


(GeV) 


hierarchy 


SPSla 


112.46 


112.45 


h6b 


SPSla' 


113.92 


113.92 


h6b 


SPSlb 


116.49 


116.44 


h3 


SPS2 


116.67 


116.61 


h5 


SPSS 


115.59 


115.52 


h3 


SPS4 


114.82 


114.81 


h6 


SPSS 


112.02 


111.92 


h3 


SPS7 


113.04 


113.04 


h3 


SPSS 


115.03 


115.02 


h3 



Table 3: Comparison of the best three- loop prediction to the one where some expansion 
parameters are cut (see text). The last column shows the chosen hierarchy. 



scenario 


1 loop 


2 loops 


3 loops 




141.11 


133.96 


134.38 


no-mixing 


121.33 


119.72 


121.07 


gluophobic 


118.52 


117.47 


117.29 


small aeflf 


120.20 


118.33 


119.10 



Table 4: Results for M^ (in GeV) for the benchmark scenarios defined in Ref. 
tan/3 = 10 and M^ = 1500 GeV has been chosen. 



where 



For definiteness we show in Tab. [3] the three-loop prediction of M^ for the SPS points. 
We show the best prediction (M^ ) and, for comparison, the result obtained by cutting 
parameters as for Fig. M {Mf^ '^^ ). Furthermore, we indicate the optimal hierarchy as 
chosen by H3m. In all cases the uncertainty due to our approximation can be estimated 
to be below 100 MeV. 

As a last phenomenological application we consider the benchmark points identified in 
Ref. [52] in order to perform the MSSM Higgs boson search at hadron colliders. In 
Tab. mthe one-, two- and three-loop predictions for Mh are shown for the four scenarios 
ii^max)) ^ "no-mixing" , "gluophobic" and "small Oeff" . Whereas for "no-mixing" significant 
three-loop effects are observed more moderate, however, still important contributions are 
obtained for the remaining three scenarios. In Tab. |l]tan/3 = 10 and Ma = 1500 GeV 
has been chosen, however, a similar conclusion can be drawn for other values. 

Let us at the end of this Section estimate the theoretical uncertainty on the prediction 
of Mh after including the three-loop corrections. We divide the uncertainty into two 
parts: (i) the theoretical error due to missing higher order corrections and other as of yet 
uncalculated corrections, and {ii) the parametric uncertainty, mainly due to the top quark 
mass, as and the supersymmetric masses. At two-loop level, a thorough investigation 
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of the theoretical uncertainties (i) has been performed in Ref. [8] (for earher work, see 
Ref. [H]). It was found that missing two- loop corrections (most importantly electro-weak 
and finite-momentum effects) can be assumed to be well below 1 GeV, while the as of 
then unknown three-loop effects could amount to 2-3 GeV, as indicated by the variation 
of the two-loop results with the renormalization scale. 

With the ata"^ corrections of our calculation, we should therefore be able to reduce the 
theory uncertainty. Instead of renormalization scale variations, however, we want to adopt 
a more conservative attitude and assume a geometric progression of the perturbative series. 
As can be confirmed at two-loop level, a conservative estimate of the uncertainty is 50% of 
the difference to the next lower order. Thus we assign 50% of the three-loop contribution 
to M/i as a theoretical error. For the MSUGRA scenarios this leads to an uncertainty 
of about 100 - 200 MeV for mi/2 = 100 GeV and to about 1 GeV for mi/2 = 1 TeV. 
These estimates cover also corrections from renormalization group improvements which 
at two-loop order lead to shifts in M^ of the order of a few GeV [?] . 

Earlier in this section, we have identified two more contributions to the uncertainties of 
type (i): The corrections beyond the quartic top quark mass contribution, A'''^™M^ , 
has been estimated to about 100 MeV and the uncertainty due to our approximation 
procedure also amounts to at most 100 MeV. Both contributions are smaller than the one 
due to missing higher order corrections discussed above. 

The parametric uncertainties can be easily estimated with the help of H3m. For definiteness 
we adopt in the following SPS2 and vary Mt and as{Mz) as follows 

Mt = 173.1 ±1.3 GeV, (22) 

as{Mz) = 0.1184 ±0.0020. (23) 

In the case of the top quark mass we observe for mi/2 = 100 GeV a variation of Mh by 
about 350 MeV which increases to SMh = 1 GeV for mi/2 = 1 TeV. The corresponding 
numbers for the uncertainty in as{Mz) read 80 MeV and 600 MeV. A further uncertainty 
is connected to the (unknown) supersymmetric parameters which may also be of the 
order of a few hundred MeV. Assuming, e.g., an uncertainty of 10% for m^^ in the range 
between 200 and 800 GeV and adopting the remaining parameters from Eq. (12T1) leads to 
an uncertainty of at most 500 MeV for Mh. Note that the parametric uncertainty is of the 
same order of magnitude as the theoretical error due to missing higher order corrections 
as estimated in the previous paragraph. 



6 Conclusions 

In this paper we have discussed the lightest Higgs boson mass of the MSSM to three-loop 
accuracy. At this order an exact calculation is out of range and thus we heavily exploit 
the methods of asymptotic expansion in order to provide precise approximations. This 
procedure has been successfully tested at two-loop order against the full result. The result 
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is expressed in terms of DR parameters for the quark and squark masses for which we 
have found very good convergence of the perturbative expansion. We provide a user- 
friendly Mathematica program H3m which allows the computation of Mh in a simple way. 
In particular, it is possible to apply various SUSY breaking scanarios, invoke a spectrum 
generator, and use the output in order to compute Mh- 

As already mentioned, with the help of the asymptotic expansion we have implemented 
analytical results valid for various hierarchies in the supersymmetric masses. H3m is set 
up in such a way that it is straightforward to include further hierarchies in case they are 
needed for future investigations. 

We have performed several studies with the new three-loop corrections. In particular, we 
have considered their renormalization scheme and scale dependence, and their numerical 
effect in the various SPS scenarios. Furthermore, assuming MSUGRA, we have considered 
the mo-mi/2 parameter space and computed two and three-loop corrections. On the basis 
of these investigations we estimate the remaining theory uncertainty of Mh to about 
200 MeV for mi/2 = 100 GeV and to about 1 GeV for mi/2 = 1 TeV. 
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Appendix A: DR counterterms 

In this appendix we provide details about the computation of the counterterms for the 
gluino, top quark and squark masses and the mixing angle in the top squark system. 

To define our framework for the computation of the top squark masses and mixing angle 
counterterms, we start from the bare Lagrangian containing the kinetic energy and mass 
terms 

cf = '-dMArd^[i^) - i^{ii:r^r\M]r [^f^^ , (24) 

where the superscript (0) labels the bare quantities, t^ and Ir denote the interaction 
eigenstates and the top squark mass matrix was defined in Eq. (JHl). 
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The top squark mass eigenstates are related to the interaction eigenstates by 

r \(0) / r x(0) 

where the matrix 7?.^ is defined in Eq. ( fTOl) . 

The wave function renormahzation can be written in matrix form 

where we have Z~^^ = l + 0{as). Thus, Z^f = l + 0(a,), Z^i"^ = l + 0(a,), Z^f = C(a,) 

and Z21 = 0{as)- 

Similarly, the renormalized mass matrix can be parametrized as follows 

H ^\=( "^ll^mii "^12^mi2 \ =M, (27) 

where the mjj stand for renormalized mass parameters. 

In order to extract the renormahzation constants it is convenient to consider the renor- 
malized inverse top squark propagator which reads 

tS-\p') = V\Z\'')'Z\I' - {Z\I')\M - nv^y^z]!' (28) 

where S(p^) denotes the matrix of self energies constructed by the ti and ^2 fields. 

To determine the renormahzation constants introduced above one has to specify the renor- 
mahzation conditions. In the DR scheme the finite parts of the renormahzation constants 
are set to zero and the coefficients of the e poles are obtained from the requirements 






= 0, 
pp 

2 



where "pp" stands for the "pole part". The conditions in the second line of Eq. 
ensure that the renormalized fields ti and ^2 are the mass eigenstates of the renormalized 
mass matrix. Accordingly, we can identify ^mn and Zm22 with the DR renormahzation 
constants Z^- and Z„^. defined through 

(m? )(°) = m|Z„.^ , z = l,2. (30) 

In the DR scheme the wave function renormahzation constants are independent of the 
masses, so we can evaluate them solving Eqs. ( p9|) for mi^ = m^^ = 0. In addition, we can 

1/2 

exploit the symmetry of the matrix Z^ and reparametrize it in terms of the common 
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wave-function renormalization constant of the fields ti and ^2 and the counterterm of the 
mixing angle, as follows 



,1/2 



z\" ^ z^ 



1/2 



COS bQt sin bQt 
- sin bQt cos bQt 



(31) 



We solve Eqs. (12^ iteratively by inserting the perturbative expansion of all quantities, 



Z = l + ^5Z«+f^V<5Z(2) + 0(a3 



vr 



-M 



vr 



a. 



vr 



vr 
aA2„(2) 



6ei'^ + 0{al) 



^sg)+(^)-Sg^ + (9(ar 



vr 



^,i = i,2 



(32) 



and take into account the reparametrization of Eq. (13T]) . At one-loop order we get 



'^^-mlSSZ^ + SZ^)+p^6Z^ 



pp 



pp 



0, 
0. 



z = l,2. 



In a first step we solve the first equations with mi^ = mi^ = to obtain 



5Z. 



(1) 



^(1) 



pz 



^'^ 



1 or i = 2 . 



(33) 



(34) 



pp 



^(1) 



Afterwards, we determine 5Zml such that the equations are fulfilled also for finite tti^.. 
The renormalization conditions for S12 and S21 in Eq. ( l29l) provide an expression for 



50^': 



(1) 
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Similarly, at two-loop order we have 



^(1) 

"12 



2 2 

rrij — mf 



(35) 



pp 



(2)^v(1)a^(1)_^2 



slf + sl^^^z, 



|(<5Z<^) + 5Z«5ZW+5Z^)) 



+(-1) 



i+lA/3(l)[ 



^(1) 



3(1)^^2 



2XV(2) 



[-2Sl^2^+(JC^(m|^-m4)]+p^5Z^ 



1,2) 



pp 



^(2) 



^(l)rv(l) 



(l)/v(l) 



sr2^ + Ei^,^5Zr+5C^(Sri^-S 



^(1) 
22 






3(2) 



(m|^-m4)(5C^ + 5C^5Z 



(i)xv(i)^ 

2 y 



(36) 



pp 
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In the following we list the results for the squark mass counterterms and the mixing angle, 
however, we refrain from providing explicit results for SZ2 and 6Z2 , as we do not need 
them for the actual computation. 

Our results for the renormalization constant of the top squark mass ttij-^ read 

1 



m 



{SZ^^ = Cp f -m? - ml + m-gmtSt + 



?Ti? — 777,? 

4 * / e 



-1^^, 



cl 



C^mjml ^ {l + C^KK - <) + 8^? (1 + C^)m,m,S, 



1 1 

vfij — vn- 



16 



+ CaCf 



9m? S^A^l - m|) + Ami Sm^mtSt 



+ CpNfTF 



32 

2 



-3m? 52(m? -m?J+4m2 m^mt^t 



+ 



+ {cl 



3m? ^?(m? -ml)+Am^ 



5 _|_ * ^ ti 



16 2 

t^y ' --"t m-gmtSt 



16 

2/^2 



+ CaCf 



-11m? 52(m|-m|)+4m2 3^.^^^^ 



32 

,2 



4 



3m? 5/ (m?^ - m| ) + 8m| + 4m,^ m^m^^^ 



16 



3m? ^2(^1 - m|) + 4m| + 4m| - 4m? m-gmtSt 



+ 



16 



+ Mt[-CACF^- + CFNfTF-^-^ 



(37) 



where we have introduced the abbreviations A'^^ = Nt + Ng, St = sin26'( and Ct = 
cos26t. Nq denotes the number of light quark flavours and takes in our case the value 
Ng = 5. Nf = 1 has been introduced for convenience. Furthermore M^ denotes the on- 
shell renormalized e-scalar mass. The corresponding results for m^^ can be derived from 
Eq. fl37p by interchanging ttt,^-^ and ttt,^-^ and changing the sign of 6t ■ 

We have employed two different approaches for the renormalization of the e-scalar mass. 
In one approach we renormalize it on-shell and choose M^ 7^ 0. In order to decouple the 
unphysical parameter M^ from the physical observables, i.e. M^ in our case, we modify 
the top squark masses by a finite counterterm. This renormalization scheme is equivalent 
to the DR' scheme [531IM] . However, in the original version of the DR' scheme the e-scalar 
mass is renormalized minimally, whereas we choose the on-shell scheme. Thus, the finite 
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counterterms we found differ starting from two-loops from the ones given in Ref. [5l] by 



Snir 



m: 



2 /^^^ 

TT 



+ CfN.Tf- 



— + CaCf ^ 



^^CpNtTp- ^^^'^^ 



4 



1,2, 



(3^ 



with L^a; = ln(/i^/?n^) and L^^ = ln(/i^/Mg^). In the second approach we set the on-shell 
e-scalar mass to zero. Of course, no finite counterterm is needed in this case. The resuhs 
for Mh obtained employing the two methods are identical, however, the second approach 
is computationally less involved. Thus for the practical calculation we set the on-shell 
e-scalar mass to zero. 

Let us mention that strictly speaking our renormalization coincides neither DR nor DR' 
due to the different treatment of the e-scalar mass. Nevertheless we use the nomenclature 
"DR scheme" . 

Finally, for the mixing angle we have 



( 2 2 N e/)(i) n n I ^tym^ 

(m^-^ - m^JdUl = CpCt \ rrigmt - 



t2' 



(ml-ml)6el 



+ CfCaCi 



(2) 



ClCt 



^F^t 



is^ - c? 



4 

rrigmt 



Stiml-ml] 



16 



-Snigmt 



3Sf(m'^ — m 



+ ClCt 



4 
rrigmt 



+ 



ti 



2 ^^ 



32 

St{ml-ml] 



+ CFNfTFQ 



rrigmt 



16 



Stimj — m? 



+ CFNjTFCt 
CFCACt 



rrigmt 



2Stmjmt ' 
m| — m? 

Stimf — -m? 



t2' 



2 16 

3mgmt 3gf(m| -m|) 
4 32 



t2' 



16 



(39) 



The two-loop counterterms given in Eqs. ( l37j) and ( l39l) can also be derived from the more 
general results available in the literature [551 - I57] . 



The on-shell renormalization constant of the e-scalar mass to one-loop order is given by 

(M2)(o) = M^ZMi , Zmi = 1 + -SZ^^ + 0{al) , (40) 
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where 



M!SZ% 



"^ml + 2ml) - ^(Me + 2ml) ' ^(^^ + < + m| - 2ml)] \ 



Ca 

^-[{2 + L^,)Ml + 2{l + L^^)m^ 



+ -^[{Q + AL^,-L^g)Ml + 2{l + L,rn)m 
N,Tf 



LV" ' —/it "Atyy-e ' -v^ ' —lJ.gl"''g\ 



2 

N T 
\^[L,tM^, + (1 + VJm| + (1 + L^,-Jm| - 2(1 + L,,)ml] . (41) 

The two-loop renormalization constant for the top quark mass in the DR scheme is known 
for long time [58l[59]. For the convenience of the reader we quote the results which are 
given by 



mj 



(0) 



m,Z„, , with Z,,, = 1 + ^5Z« + (^) 5Z^^] + 0{al) , (42) 

and 



(0) 



with 



*zS = -c,i. 




^ZS - (icHf.C.C.-lc.A^.T.)^^ 




\o io o ye 


(43) 


the one-loop eluino mass counterterm 581159] defined through 




= m^Z™-, with Z^_ = 1 + ^<5Z« + 0(a2) , 


(44) 


^z^. - [-\^^^\^fT:)\. 


(45) 



Appendix B: Modification of the DR scheme: MDR 

As discussed in Section [31 the renormalization constants of the top squarks have mass 
dimension two. Thus, for some hierarchies considered in this paper parametrically (and 
numerically) large corrections might appear which are absent in the on-shell scheme. In 
order to cure this problem we introduce additional finite corrections in the corresponding 
renormalization constants which ensure that the renormalized result for the Higgs boson 
mass is free of these potentially dangerous contributions. 
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In the following we provide analytic expressions for the finite shifts introduced in the top 
squark mass counterterms as compared to the DR scheme. According to the discussion 
in Section [HI one can distinguish three cases for the mass hierarchies. 



Case (i): 



mg >mi-^, 



mi. 



1,2) 



a^\2 



IT 



CfN^Tf 



2 

mf 



L,^ + C(2) 



(46) 



The label Ng = 5 has been introduced for convenience and for the logarithms the abbre- 
viation L^g = ln(/i^/m?) has been introduced. 



Case fiil 



mi2 > ^i 



1 1 or 
tl 

mz 



1 m| 

F^F — 5- 



CpTp 



1 

'4 



2 -^Mt2 



:C(2) 



(47) 



In this equation we have L^^-^ = ln(/i^/m 



t2' 



Case (iii) 



m-g -:^mi^, 



1, 2) and rriq ^ rrig 



m 



MDR' 



ntf 
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n mi \ n 
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21 
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-g - -^L^g - -L^g - -L^gLf^g + " L^ - + 2^(2) J ■^;;^ 
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43 5 
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4 
^^ 
2 2 

m^rrij 



67 



288 



1l- 

24 ^^ 



?7l? 



4 2 
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(48) 



Here A^, 



1, -^/xg 



ln(/i^/m?) and Lgg = ln(m|/m?). 
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Case (iv): rrig ^ tti^. , (i = 1, 2) and rriq ^ nig 
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MDR' 
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rar 



1 + "^Cf [1 + L,g] ^ 
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TX 



cl 
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21 7 9 2 
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2 H 



(49) 



All the masses on the r.h.s. of the Eqs. ( I46ll . ( I47ll . ( I48ll and ( l49ll are DR masses. Cases (i) 
and (ii) are applied for all hierarchies whereas cases (iii) and (iv) are only used for (h6) 
and (h6b), respectively. Let us also mention that the above formulae are valid for the 
case Mf = 0. 



Appendix C: SPS scenarios 

In the following table we list the input values for the MSUGRA SPS scenarios as defined 
in Refs. [19l[20]. All masses are given in GeV and sign(/isusY) = 1- 







Points 




Slopes 




mo 


mi/2 


Ao 


tan/3 


niQ 


^0 


SPSla' 


70 


250 


-300 


10 


— 


— 


SPSla 


100 


250 


-100 


10 


OAmi/2 


-OAmi/2 


SPSlb 


200 


400 





30 


— 


— 


SPS2 


1450 


300 





10 


2mi/2 + 850 





SPS3 


90 


400 





10 


0.25mi/2 - 10 





SPS4 


400 


300 





50 


— 


— 


SPSS 


150 


300 


-1000 


5 


— 


— 
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